Substrate gating of contact resistance in graphene transistors 
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Metal contacts have been identified to be a key technological bottleneck for the realization of 
viable graphene electronics. Recently, it was observed that for structures that possess both a top 
and a bottom gate, the electron-hole conductance asymmetry can be modulated by the bottom gate. 
In this letter, we explain this observation by postulating the presence of an effective thin interfacial 
dielectric layer between the metal contact and the underlying graphene. Electrical results from 
quantum transport calculations accounting for this modified electrostatics corroborate well with the 
experimentally measured contact resistances. Our study indicates that the engineering of metal- 
graphene interface is a crucial step towards reducing the contact resistance for high performance 
graphene transistors. 



Since its experimental isolation [IH3] , graphene has at- 
tracted significant attention from the scientific commu- 
nity, due to its unique electronic structure and physical 
properties [U [5] . Its excellent transport properties and 
the ability to tune the carrier concentration with elec- 
trical gates, also makes it a material with great techno- 
logical promise. Potential applications range from RF 
devices and transistors (6j-(8 to bio-sensors ^ and flexi- 
ble electronics [TU]. The metal- graphene contact, is how- 
ever, a key technological challenge for graphene-based 
electronic devices. For current-generation silicon metal- 
oxide-semiconductor field-effect transistors (MOSFETs), 
the International Technology Roadmap for Semiconduc- 
tors calls for a resistance of 80 Vl-fim per contact, which 
is about 10% of the transistor's on-resistance Vdd/Iqn 
[llj . Graphene's excellent transport properties should 
produce transistor on-resistances considerably lower than 
those of silicon MOSFETs. To realize the performance 
potential afforded by the excellent transport properties 
of graphene, exceptionally low contact resistances will be 
required [SI IZ]- It is, therefore, essential to develop a 
thorough understanding of metal-graphene contacts and 
of the fundamental lower limits for the contact resistance. 
In this paper, we develop a model that explains the re- 
cently observed substrate modulation of contact resis- 
tance in graphene transistors [12l US] . We argue that this 
effect is due to the presence of an effective thin metal- 
graphene interfacial dielectric layer. Using this model, 
we estimate two important components of the series re- 
sistance and establish lower bounds for the contact re- 
sistance. The study provides an improved understanding 
of the metal-graphene contact that may prove useful for 
improving device performance. 

Graphene is sensitive to external perturbations due 
to its all-surface and zero volume nature [5]. Charge 
transfer between metal and graphene due to a workfunc- 
tion difference dopes the underlying graphene |14j . Con- 
tacts, therefore, introduce a built-in electrostatic junc- 
tion within graphene, which was observed experimentally 
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using scanning tunneling microscopy jl5| . A distinct ex- 
perimental signature was the asymmetry of resistance in 
back-gated devices [16l [17] . The sign of this asymmetry 
reflects the doping of the graphene underneath the metal. 
Recently, devices with top and bottom gating schemes, 
as shown in Fig.jl] were experimentally realized [T^ 113) . 
Measurement of resistance vs. top gate voltage (Vr) re- 
vealed an asymmetry, whose sign and magnitude were 
modulated by the back gate voltage (Vb). This observa- 
tion strongly suggests that the graphene doping under- 
neath the metal was substantially modulated by the back 
gate voltage. 
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FIG. 1: (a) Top/bottom gated graphene structure, and (b) mod- 
eled potential profile across metal coated part and channel part of 
graphene sheet. With Ct/(Cb) we denote the top/(bottom) gate 
capacitances, whereas Cm accounts for the interfacial capacitance 
between metal and graphene. The positive/(negative) sign of the 
difference Bp — E]j between the Fermi level Ep and the Dirac 
surface Eo accounts for n/(p) doping of the respective region. 
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It is instructive to recall that in standard metal- 
semiconductor junction theory, the Schottky barrier 
height is given by the difference between the metal 
workfunction, 4>m, and the electron affinity, xs, of the 
semiconductor, when there is no Fermi level pinning. 
However, the metal-graphene binding is comparatively 
weaker. The nature of the interfacial metal-graphene 
chemical bonding is still a subject of theoretical study 
[IB] . Electrically, this interfacial layer can be modelled 
with an interfacial capacitance Cm- If Cm ^ Cb and 
Cm > Cqm ~ e'^DlEp — E^i), where D stands for 
the density-of-states, then one expects the Dirac point of 
graphene Ejji to be stationary with respect to (pM- The 
above limit describes the standard metal-semiconductor 
junction, since Cm is generally significantly larger than 
the semiconductor capacitance. HCm < Cqm, it is possi- 
ble to modulate Eui with an applied voltage Vb- Consid- 
eration of this fact is key to explaining the experimental 
results [I2J of the device shown in Fig.[l] 

The model: From Gauss's law, the electrostatic equa- 
tions governing the graphene underneath the contact (re- 
gion 1) and top gate (region 2) are 

Cm{Edi + H) + Cb{Edi - eVB) = e'n^ (1) 
CT(Eu2-yT)^CB{Em-eYB)^e:'n2, (2) 

where b(\) = 0m — <^q is the workfunction difference 
between metal and graphene, and eui — E\j^li:(hvf)^ . 
Note that Ep is taken to be zero as reference. The tran- 
sition between region 1 and 2 is described by an ana- 
lytical screening model [19] . assuming a linear graded 
junction. The junction resistance, TZjunc, is then com- 
puted quantum mechanically using a previously devel- 
oped mode space non-equilibrium Green function method 
for graphene ^D]. For our ballistic study, the effect of 
temperature impacts our results only through the ther- 
mal smearing due to Fermi-Dirac distribution function, 
which is included in our study. 

In the experiment |12| . the gate capacitances are 
known, i.e. Cb ~ 1.15 x 10"'' Fm'^ and Ct = 42Ci3. 
On the other hand. Cm is a quantity to be determined. 
Ti/Pd/Au is used for the metal contacts. For our calcu- 
lations, we assumed Sip « 25 meV, a value that is sensible 
for our experimental metal stack [37]. Fig. [2^ plots Edi 
as a function of Vb and Cm- As expected, Vb modulates 
the doping in region 1, with greater ease when Cm is 
smaller. In the experiments, it was observed that the con- 
ductance as function of Vt exhibits the least asymmetry 
when Vb ~ lOV (see Fig. [3^). This suggests that when 
Vb ~ 10 V, Ebi ~ (the charge neutrality point), which 
then also allows us to pin down Cm to he Cm ~ 400Cb. 
On the other hand, a different choice of 5<j) would corre- 
spondingly yield a different Cm, as illustrated in Fig.jSja. 

Results: The measured conductance, G{Vt), for dif- 
ferent values of Vb is shown in Fig. [3^. The observed 
asymmetry in C{Vt) changes sign at about Vb = lOV . 
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FIG. 2: (a) Surface potential (left), and (b) modulation of the 
neutrality condition (right) of the metal coated side of graphene. 
The doping of the metal coated side of graphene is modulated with 
the backgatc Vb and depends also on the value of the interfacial 
capacitance Cm and the work function difference 5<j). The crossing 
line(s) in both figures refiect(s) zero surface potential (Dirac lines). 
In the inset (right), the modulation of the neutrality condition of 
the metal coated side of graphene is applied, in direct correspon- 
dence to the equation Cm^4' ~ CbVb = 



To facilitate comparison between experiment and a bal- 
listic theory, we extract the odd component of the resis- 
tance [3T] from experiment, given by TZodd = \ V^{5Vt) — 
TZ{—5Vt)], where 5Vt is Vt with respect to the Dirac 
point voltage. The quantity TZodd then allows for quan- 
titative comparison between the experimentally mea- 
sured G{Vt) and the numerically calculated TZjunciVx)- 
Although interface charge, moisture, and chemicals in 
the vicinity may impact the transport properties of the 
graphene devices, these effects can be minimized by care- 
ful control of the fabrication process and the measure- 
ment conditions so that any asymmetry can be attributed 
entirely to the graphene p-n junction. 

In Fig.[3|D, the striking experimental observation of 
asymmetry inversion, as it is observed in TZodd vs Vb = 
—AOV and Vb = 4:0V, is compared with our modeled 
TZodd for a work function difference 6(j) = 0.025 eV and 
an interfacial capacitance Cm = 400Cs. As elucidated 
previously, modulation of the doping of graphene un- 
derneath the metal is possible because Cm is not large 
enough to completely dominate over Cqm- In fact, for 
a moderate carrier concentration of 1 x 10^^ cm~^, one 
obtains Cqm « 0.75Cm- If one assumes that the metal- 
graphene interfacial layer is an air gap. Cm = 400Cb 
would then translate to a physical thickness of only 2 A. 
This is in good agreement with recent density functional 
studies j22j , with a predicted metal-graphene binding dis- 
tance of ~ 3.5 A. The possibility of sign inversion is, 
however, conditional. For example, if we use a value of 
Cm = 4000Cm instead of Cm = 400Cb, sign inversion 
of TZodd would not be observed within the Vb range of 
interest, as shown in Fig. [3]:. In Fig.|3ji, we show that 
the increasing odd resistance with increasing \Vb\, as it 
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FIG. 3: Study of the conductance asymmetry of the Chen/Appenzeller experiments |12| . (a) Experimental G — Vt curves for different 
back gate voltages Vb, ranging from Vb = —40 V to Vb = 40 V. (b) Comparison of model prediction versus experiment for the odd part 
of resistance {"R-odd) vs top gate voltage Vt for Vb = —40 V and Vb = 40^, assuming an interfacial capacitance Cm = 400Cs, and a 
temperature of T = 300 -fC. (c) Model prediction of (liodd) for the same voltages for an elevated interfacial capacitance Cm = 4000Cs, 
for T = 3Q0K (solid line) and T = OK (dashed line), (d) Odd part of resistance (Tlodd) vs Vt for Vb = OV and Vb = -40 V with 
Cm = 400Cs at T = 300 K. (e) Experimental R — Vt curves for Vb = —40 V and Vb = 40 V, versus our NEGF model, assuming an 
interfacial capacitance Cm = 400Cfl, a work function difference 50 = 0.025 eV and a temperature of T = 300 -ft'. 



is observed in the experiment for Vb = 0, — 40 can also 
be captured with our simulation, using Cm = 400Cb. 

Up till now, we only considered TZjunc and its con- 
tribution to the asymmetric part of the contact resis- 
tance. Previous studies [23l El] of the intrinsic transport 
properties of graphene on substrate allow one to make 
reasonable estimates of the channel resistance, TZcha, by 
including contributions due to acoustic/optical phonons 
and substrate-induced remote phonons. However, their 
contributions relative to TZjunc sn'e not as significant 
|38j . A modest mobility of ~ 500 crn^ /Vs was extracted 
from the experiment in vicinity of the Dirac point, sug- 
gesting high levels of impurities |25j. Here, we model 
the impurity limited resistivity with the Laudauer for- 
mula using a mean-free-path proportional in energy i.e. 
p^i = h/2e^ X M-'^[L-^ + {aE)-'^], where a is used to fit 
the mobility and M is the number of modes normalized 
to W. In Fig. [3^, we compare the experimental resis- 
tance with the calculated sum TZc = Tlcha + TZjunc- By 
construction, our model does not capture the physics at 
the Dirac point. However, far away from the Dirac point, 
one observes an unaccounted for excess resistance in the 
experiment of ~ 500 the origin of which is the subject 



of subsequent discussion. 

Other contributions to contact resistance include the 
current crowding effects due to access geometry and the 
presence of the metal-graphene interfacial layer. The 
latter implies that current has to tunnel across a di- 
electric layer, encapsulated in the electrical quantity 
Cm = 400Cb. The tunneling resistivity ptun can be 
estimated using a model for quasi-bound electrons |26) . 
commonly used in the study of gate leakage current 
in semiconductor inversion layers |27| . It is given by 
Ptun = Pr/e^D, where P is the Wentzel Kramers Bril- 
louin tunneling probability estimated to be « 0.3, and 
T is known as the classical bounce time [26], which in 
graphene is simply r = 2tg/vf, with tg and Vf being 
the thickness of graphene and Fermi velocity respectively. 
Here, we take tg to be the carbon-carbon bond length, 
1.44 A. This yields us a tunneling resistance estimate of 
Ptun ~ 5.2 X 10"^ il-cm? for a pair of contacts. This value 
is reasonably close to a recently reported [3S] experimen- 
tal value of « 5 X 10"^ Vt-cm?. 

Expanding the prior analysis to the spatial distribu- 
tion of carrier flow, the current crowding effects due to 
the metal-graphene acesss geometry lead to an effective 
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FIG. 4: Estimates of various components of contact resistance in 
graphene transistors. With the long-dashed line, we denote our 
modeled ballistic prediction for the junction resistance plus the 
tunneling resistance TZc = 'R-junc + "R-tun, for a pair of interfaces. 
After subtraction of Titun, we get the ballistic result for the lower 
bound of the contact resistance TZjunc (short-dashed line), assum- 
ing a work function difference (f>M — <I>G = 0.1 eV. Last, an in- 
crease of the work function difference to (I>m — <j>Q = 0.3 eV reduces 
further our ballistic predicition for the junction resistance TZjunc 
(solid line). (The region near = should be disregarded since it 
is strongly affected by spatial potential fluctuations, which are not 
considered in our model.) 



electrical area for the contact where current flows, given 
simply by X Lt, where Lt = a/ Pc/ Pg is commonly 
known as the transfer length. With pc and pg , we denote 
the specific contact resistivity and the sheet resistance 
of the graphene layer underneath, respectively. Using 
our estimated ptun ~ 5.2 x 10"^ ^-cw? and a sheet re- 
sistance Pg = 1660 f2/n from the experiment [12], at a 
channel carrier concentration ris ~ 5 x 10^^ crrf^, we find 
Lt ~ 560 nm. Given the experiment to experiment vari- 
ations in Pc and pg , this value is within a reasonable range 
[2H]. 

Perspective: Recently reported contact resistance 
values lie in the range 7^c « 600 - 10^ n-pm [H [21- 
l32] (see Suppl. Info.). These values are considerably 
above what is required for high performance transistors 
[llj . In this section, we examine several issues related to 
the fundamental limit to TZc- 

In Fig.|4l we consider the contact resistance for a pair 
of interfaces TZc, incorporating both the tunneling com- 
ponent TZtun ~ 520 that we extracted previously, and 
the junction component TZjunc assuming a work function 
difference (f>M — (j)G = 0-leV (long-dashed line). It is 
expected that TZtun can vary considerably from experi- 
ment to experiment, due to different contact materials, 
interface conditions, etc. [T51 [^5H55] . Significant reduc- 



tions in TZtun will be crucial, and the contact metal and 
deposition and annealing conditions will be critical fac- 
tors. Recently, it was demonstrated that a low power 
plasma O2 treatment prior to metal deposition is bene- 
ficial in improving TZtun ^34 . Another very recent study 
demostrates that TZtun is temperature dependent, possi- 
bly due to an expansion of graphene-metal distance |33j . 

The ballistic junction component TZjunc, on the other 
hand, seems to have a more universal nature, since it is 
limited by the electrostatics condition, and the number 
of conducting channels bottleneck [20 . To illustrate this, 
we plotted TZjunc for (fiM 4>g — 0.1 and 0.3 eT^ in 
Fig.|4] For channel carrier concentration < 0, a lower 
resistance plateau forms. The resistance value of this 
plateau is described by the quantum contact resistance 
Rq = h/2e^ X M~^W [35], where M is the number of 
current-carrying modes, and this limit is imposed by A0 
at the metal side of the junction. Indeed for A0 = 0.1 eV 
(short-dashed line), one finds Rq « 12>AQ,-pm, while for 
= 0.3 eV^ (solid line) one finds Rq ~ AhVt-pm. In 
the latter case, however, for moderate negative values of 
Us, where the plateau hasn't been reached, quantum re- 
sistance Rq is limited by the number of modes in the 
channel, and Rq = h/2e^ x l/2-\/7r7n7 holds valid in- 
stead. The elevated right branch, when Ug > 0, is due to 
the effect of interband tunneling. Selection of appropri- 
ate metal workfunction or approaches to engineer sharp 
p-n junction, i.e. chemical doping [5], are promising di- 
rections. 

In conclusion, we have proposed a model that ex- 
plains the gate dependent resistance asymmetry observed 
in experiments, and provides increased understanding 
of the different components of the contact resistance in 
graphene transistors. We show that the existence of an 
interfacial layer between the metal and graphene can ex- 
plain the back gating of the contacts observed by Chen 
and Appenzeller jl2j . The importance of such an in- 
terfacial layer to devices was also pointed out in two 
very recent studies. Robinson et al. show that resid- 
ual photoresist following lithography can lead to high 
contact resistance [34], and Xia et al. explain the tem- 
perature dependence of the contact resistance in terms 
of the graphene to metal distance [33] . Our work and 
these studies point to the need for further understanding 
of the metal-graphene interfacial layer as a prerequisite 
for engineering a good contact resistance for graphene 
electronics. 
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